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The Fermat principle indicates that light chooses the temporally shortest patri. The action for 
this "motion" is the observed time, and it has no Lorentz invariance. In this paper we show how 
this action can be obtained from relativistic action, and how the classical wave equation of light can 
be obtained from this action. 
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I. MASSLESS LÍMIT OF RELATIVISTIC 
PARTICLE AND DIMENSIONAL REDUCTION 

Let us start with the action of relativistic massive par- 
ticle, i.e. the proper time 



S'ò 



m / dr^J x^x l 



(D 



But we can not take massless limit of this action. So we 
uso the diffcrcnt action which is equivalent to S'ò 



Si 



1 



dTyfg(g 1 x fl x fl + m 2 ), 



(2) 



where g is the auxiliary field. Let us show the equivalence 
of these two actions. The Euler-Lagrange equation for g 
is 



.9 



m- 



(3) 



By putting above equation for g into the action Si, we 
obtain So- So we observe the equivalence. To consider 
massless limit, we use Si instead of So- The action for 
massless particle is then, 



Sa = — - / drg 



-1/2^, 



Two Euler-Lagrange equations are 



dT y/9 



x^if, = 0. 



The latter equation gives 

x° = ±V¥. 
Then the first equation of 10 becomes 



(4) 



(5) 



(6) 



The first equation shows that ratio of g and x 2 is con- 
stant. And so by using constant A, we can write 

5 = Ax 2 . (8) 

By putting this into the second equation of J7J, we have 



d x l 



dT 



0. 



(9) 



The action which gives © as Euler-Lagrange equation is 
S 3 = I J drVÏ^, (10) 

where the speed of light c is introduced for the action to 
have the dimension of time. S3 is the time along the path 
of light, and it is really the action expressing the Fermat 
principle. Though it seems non-relativistic, it should in- 
troduce the wave equation of light. This point will be 
made clear in the following sections. Note that the di- 
mension of configuration space in S3 is 1 dimension less 
than the one in So- The massless limit induces the dimen- 
sional reduction in action, which is similar to the change 
of physical degree of freedom: from massive vector field 
(3 degrees of freedom) to massless vector field (2 degrees 
of freedom). Another point is that action S3 still have 
time-reparametrization invariance Due to this 

local symmetry, the method to obtain the wave equation 
becomes bit complicated. This is discussed in the next 
section. 

II. TIME REPARAMETRIZATION SYMMETRY 

Let us consider the problem to obtain the wave equa- 
tion of light from action S3. But S3 is not enough to 
express the Fermat principle. Because the minimisation 
of S3 can also be understood as minimisation of path- 
length. To make clear the meaning of minimisation of 
time along the path, we generalize the action with re- 
fractive index field as 



S oc "time" = — 
c 



7 \(x)[(§n/ 2 dr, (ii) 
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where c is the speed of light, n(x) is the refractive index 
and t is only the parameter but not time. íít[(^) 2 ] 1 / 2 is 



2 



the infinitesimal spatial length along the light path. This 
action is invariant under the transformation 

t--/(t), 

if the conditions df/dr > 0, = /(0), and T = f(T) 
hold. This is well known as time reparametrization y|,|2|. 
The equation of motion as "particle mechanics of light" 
is 



-j-[n(x) —- ] — Vi 2 Vn(i). 
dr ^/x 2 



(12) 



This equation can be re-expressed by using the length 
parameter s; ds — Vïc^dT, 



d 2 x 1 rr ± rr ± dx. dx. 
= [Vn-(Vn-— ) — ] 



ds 2 



ds ds 



(13) 



The right hand side can be interpreted as some-kind 
of "force" . The property of this force is the followings. 
When Vn./ /(dx/ds), then Vn =| Vn | dx/ds, and the 
force vanishes. When Vn _L dx/ds, then r.h.s. is Vn/n. 
Further, (|lll) can be expressed as 



S 



^jgijX l xHT, 



(14) 



with conformally flat mètric: = (n/c) 2 Sij. Then the 
equation of geodèsic line is the following. 



d 2 x l 
ds 2 



i dx j dx k _ 



(15) 



where ds = 7 -^^-ds. This equation is equivalent to (|12Jl . 
where the Christoffel Symbol is 



r =hs i — + s i — 



dn 



Our program is the following. To consider the quanti- 
zation procedure, it is necessary to construct dimension- 
less action. This is done by multiplying the frequency of 
light: lo (monochromc light) . We have 



5' 



n(x)[(f ) 2 Y' 2 dr. 



Then we consider the path- integration, Q.Q 
V[xf,3Bi,<J\* = elSlC] - 

pathC 



(16) 



(17) 



To obtain the time which we observe, we sum up the 
function í with various freqüències. 



®[x f ,t: x l ,0] = — J duj^ 



Xf,Xi,u!\ st e 



(18) 



Then we seek the wave equation which [xf , t : Xi,0] 
obeys. Note that we have no H in the path-integration. 



This is because we are treating Light: classical wave in- 
stead of Schròdinger equation. 

The explicit form of path-integral 1)17(1 is the following. 



9, 



T>x exp [í - 



i{x){Q?Y' 2 dr]. (19) 



But the treatment of this style action is not convenient. 
So we change the action into the canonical form. From 
the original action S, we can define the canonical mo- 
mentum Pí(t) which is conjugate to x l . Then we have 
the Hamiltonian constraint Ti. = p 2 — {nuj/c) 2 , and the 
canonical action. 



S = 



í drL(r) 
Jo 



dr [p l x i - N(p 2 -uj 2 n 2 /c 2 )' 



(20) 

where TV is the Lagrange multiplier field. The equation 
of motion from this action is the same as the previous 
one 1|12[) . Time reparametrization invariance exists also 
in this action. The transformation property for each field 
is 



x\t) 
N(t) 



Pi(f(r)), 
N(f(T))f(i 



(21) 
(22) 
(23) 



This symmetry is canonically induced from the first class 
constraint Ti. The canonical infinitesimal transforma- 
tion is given as 



5x l = {x\p 2 
Spi = {Pi, P 2 



2 /c 2 }e = 2 Pi e, 



(24) 



uj 2 n 2 /c 2 }e = 2u 2 nd l ns/c 2 .(25) 



The transformation for N(t) is defined for action to be 
invariant as 

SN = è. (26) 
The symmetry is as follows: 

SL(x,p,N) = ± [e(r) (p 2 + Lü 2 n 2 /c 2 )]. (27) 

The invariance of action requires e(0) = e(T) = 0, as the 
boundary condition which we will use later. 



III. PATH-INTEGRATION 

Under the transformation l|26|) . the naive measure for 
N is invariant. So we can safely consider the path- 
integral: 



Vx Vp VNe 



iS 



(28) 



Let us use the gauge condition: N = that fix the func- 
tional form of infinitesimal parameter as 



N = -> e = 0, 



(29) 



3 



which is the similar relation as Lorentz gauge in U(l) 
gauge theory. Then we can construct the Faddeev-Popov 
determinant as: 



A _1 [JV] = / VeS(N £ ) 



(30) 



where N £ means the gauge transformed field N(t) by 
time-reparametrization e, and T>e is the Haar mea- 
sure of time-reparametrization group. Then our Path- 
integration has the form 



VxVpVN A[N] S(N)e 



iS 



(31) 



A[N] can be calculated under the gauge condition N = 
with infinitesimal parameter e as follows: 

A[N] U =0 = [ I VeS(e)}- 1 = Det[^]. (32) 

Then we arrive at the final form 

í d 2 
= / VxVpVN [Det—] 8{N)é s . (33) 

By using the £ function regularization given in appendix 
A and boundary conditions for e, we calculate the deter- 
minant as 



717T s 



(34) 



By remembering TV > is used to construct the equa- 
tion of motion from action l|lfi|) . we present zero-mode 
integration for N. 

/>oo 

VN 8{N) = dN (35) 
Jo 

Then we have the path-integration. 
* sí = T J VxVp J dN 

exp [i í dripié - N{p 2 - u 2 n 2 /c 2 )}} (36) 
Jo 



where we defined as n 2 (x) = 1 — V(x)c 2 /u 2 . Let us 
rewrite this function as follows: 



dSe iuj2s/c2 ^ {x = x(S),x = x{0),S), (39) 



where VPo is the kerncl for non-relativistic Schròdinger 
equation with mass 1/2. 



(i-ïe + V 2 - V(x))* {x, x ,S) = S 3 (x - x )S(S). (40) 
dS 

Then we multiply e lüJ /c to both hand sides, and intc- 
grating by S, we have 



(-T + V 2 - V(x))^ st (x,x 0l io) = S 3 (x-x ). 
& 

Since V(x) = (1 — n 2 )uj 2 /c 2 , we obtain 



(V 2 + n 2 ^-)4- sí (a;, x , u) = S 3 (x - x ). (41) 



Further, by defining the time dependent function; 

*(ar, x , t)=^J e- l ^ st (x, x Q ,u)diü, (42) 
we have the result 

n Fi 

( y2 - )*( x > x o,t) = S 3 (x ~ x )S(t). (43) 

Schulman has done the reverse way. 3 He started 
from wave equation l|43l) . and changed to the path inte- 
gral formulation iJSSJ. Then he took the stationary point 
approximation to show the Fermat principle. This proce- 
dure is essentially different from our method. Our start- 
ing point is the Fermat principle itself, and by using the 
technique of quantization procedure of gauge theory, we 
come to this wave equation. If the medium has disper- 
sion, the integration by lo cannot be cxplicitly done. The 
wave equation for dispersive medium is given in integra- 
tion form. 



IV. WAVE EQUATION 

Performing the momentum integration, we obtain the 
another form of path-integration expressed only by con- 
figuration variables. There appears the infinite multipli- 
cation of N as the measure. But we can absorb such a 
measure into Dx by changing the time like variable r as 

dr' = Ndr, S = NT. (37) 

This is explicitly explained in appendix B. Then we ob- 
tain 

* rt = / dSe iuj2s ^ c2 / Vx exp[i / dr' {-± 2 -V{x)}], 

(38) 



V. CONCLUSION 

We have shown two facts. First the massless limit of 
action for relativistic massive particle becomes the ob- 
served time along the path of light, which realizes the 
Fermat principle. Second the quantization of the system 
with this action really induces the massless wave equa- 
tion as light. The introduction of Fermat Principle from 
wave equation is well known on many books, but in this 
paper we have shown completely reverse way by using 
the quantization technique for gauge theory. The gauge 
symmetry is the time reparametrization, and the careful 
treatment of such symmetry is important to obtain the 
wave equation. 
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1. Appendix A 

We show how the zeta-function regularization is per- 
formed. Let the operator A, and its eigenvalues as 
A n = an x , where n is the integer and n > 0. Then 
we calculate 



detA = Y[ A„ 



we find the relation that 



dÇ(s | A) 
ds 



detA = Yl 



(44) 



detA = exp[tr{logA}] = exp[- dC(s , 1 A) | s=0 ], (45) 

ds 

where the definition of ((s \ A) is as follows. 

oo 

Q{s\A)^Y. X n S - (46) 

n=l 

This can be written as 

oo 

((s\A) = J2a- S n- sx . (47) 

n=l 

Therefore, 



| s=0 =log(a)C(0)-xC'(0), (48) 



where Ç(x) is the usual zeta function defined as 

oo 

C(a0 = 5>-». (49) 

n=l 

We have the formula 

C(0)=-l/2, C'(0) = -(log27r)/2. (50) 
Then we obtain 

.m^^^^pi. (51 , 

ds ya 
By using above formula, we obtain 



(52) 



This is the relation we used in section 3. 

2. Appendix B 

After the momentum integration, (|36|) reduces to 
* sí = T j™ dN j-Dx ([[N- d / 2 ) 



x ex pï jf dr[^ + N(^n (53) 



where d is the spacial dimension. Then we remember the 
explicit form of path-integration by x, such as 



= T dN 



d M-l 



lim [f] TT dx l {k){—^—Y 



i=l fc=l 



M-l 



ex p * t 



fe=0 



(x(k + 1) - x{k)f 
ÏÏVe 



un. 



(54) 



where e = T/M. We define 

e = eN (dr' = Ndr). 

Then we obtain 

d M-l 



= T / diV lim 



M— ►oo 



n n ^( fc ) 



=i fc=i 



M-l 



1 JM ^— \ 

(tt- -) 2 exp i 2^ 



(x(fc + 1) — a?(fc)) 2 ^ -( wn ) 2 



4e 



JVÏ' 



= T J dN J Vx expi J dr'[— + ( — Y). (55) 
Let us define 5" = NT, dS = TdN, then we obtain 



* sí = y dsj-Dxexpij dr'^ + ^Y)- 
This is equivalent to equation (|38|l . 



3. Appendix C 

Here we show the BRST construction of path integra- 
tion (33). We start with the action (20) with symmetry 
under (24), (25), (26). Let us introduce the BRST trans- 
formation by the replacement: 



e(r) - C(t), 



(56) 



where C is the ghost field with odd grassmannian parity. 
Then we may have 

U 2 

6 B x i = 2 Pl C, S B Pt = 2^ r Cnd i n, 5 B N = C. (57) 
cr 

The BRST transformation of ghost field is determined by 
the nilpotency as 



5 tí C = 0. 



(58) 



But the related BRST charge can not be constructed, 
because BRST transformation for N field includes time 
derivative. To solve this difficulty, we must introduce 
the other fields: tt, and C, and supposing the Poisson 
brackets 



{N,ir} p = l, {C,C} P = 1. 



(59) 
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Then we write the BRST charge which governs the BRST 
transformation. 

Q' B = C(p 2 - {™f) + ttC, 6 b A(t) = -{Q' B , A(r)} 

(60) 

It is clear that the form of BRST charge is ambiguous 
since it includes "time" derivative of the field. To avoid 
this difHculty we need another kind of ghost field named 
F(r) and its conjugate field F, and require the equation 
of motion 

C = F. (61) 
Then we have the nilpotent BRST charge 

Q B = C(p 2 -( — ) 2 ) + ttF, {Q b ,Q b } p = 0. (62) 
c 

All the variables are 

x\N, C, F; Pi ,n,C,F. 
Then we have the BRST transformation for each field 

Lü 2 

5 B x l = 2 Pi C, S B pi = 2—Cnd t n, 6 B C = 0, (63) 

^C = -(p 2 -(— ) 2 ), S B N = F, (64) 
c 

S B F = 0, S b F = -tt, S b tt = 0. (65) 

It can be seen that the latter 4 variables (N, tt, F, F) form 
the BRST quartet ^ÏJ. Further the action needs to in- 
troduce the kinetic term for all the variables and also the 
equation C = F. The solution is easily found in the form: 

L = Pl x l - N(p 2 -( — ) 2 ) +7rN+FF + C(C-F), (66) 
c 

which is the reasonable form since 

L= Pl x l + irN + FF + CC-{Q B ,NC}. (67) 



Since the BRST transformation is constructed canoni- 
cally, the kinetic terms are invariant and the another 
term is also invariant from the nilpotency. Therefore the 
Lagrangian is totally BRST invariant and variation by 
C gives the required equation C — F. We should noticc 
that the gauge is automatically fixed as N — 0. Now we 
have the path-integration 



* = / VxV P VNVttVCVCVFVF 'eh LdT . (68) 



The C integration requires S(C — F), and so the F inte- 
gration changes the term 

FF -> FC. 

— ,2 

Then F, C integration gives the determinant Det4-s, 

and 7r integration induces S(N) . Therefore we again come 
to the functional 

r d 2 

* = / VxV P VN \Det—]5(N) 
J dr 2 

exp[z / d T { Pl x l - N(p 2 -( — ) 2 )}]. (69) 

This is the same equation as l|33f) . So we have the prop- 
agator of the massless particle as light. 
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